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The strong coupling phase diagram of the spinless fermions on the anisotropic triangular 
lattice at half- filling is presented. The geometry of inter-site Coulomb interactions rules the 
phase diagram. Unconventional charge ordered phases are detected which are the recently 
reported pinball liquid and the striped chains. Both are induced by the quantum dynamics out 
of classical disordered states and afford extremely correlated metallic states and the particular 
domain wall- type of excitations, respectively. The disorder once killed by the quantum effect 
revives at the finite temperature, which is discussed in the terms of the organic 0-ET2X. 
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Charge degrees of freedom provide particular orders 
and disorders under geometrical frustration, as has been 
discussed from the time of the Verwey transitions. 1 Many 
recent topics are related to those on the anisotropic trian- 
gular lattice(ATL). They are experimentally explored in 
the organic solids as a possible " non-magnetic" Mott in- 
sulator 2 and organic thyristor. 3 Theoretically, the phase 
diagram of the half-filled Hubbard model on ATL seems 
to suggest the particular "spin liquid" insulator, 4 which 
has been searched for from the past as a nonmagnetic 
state free from any kind of ordering. 5 

Quarter-filling is characterized by the charge order- 
ing(CO) phenomenon, as studied by the extended Hub- 
bard model(EHM) in the context of organic solids. 6 ' 7 
The interplay of CO and frustration is investigated in the 
phase diagrams based on ATL. 8,9 However, the two di- 
mensional systems under geometrical frustration requires 
a delicate and difficult numerical treatment, thus a clear 
answer is still not obtained. In the present letter, we sim- 
plify the problem and deal with the spinless fermions at 
half-filling to extract the intrinsic character of the charge 
degrees of freedom. We obtain a convincing phase di- 
agram which retains a strong coupling nature even in 
a rather weak coupling region. In addition to the pin- 
ball liquid reported in Ref.[10], another unconventional 
strongly correlated CO state with fluctuating stripes is 
found in the finite size cluster. Its transport character is 
sensitive to the geometry of quantum dynamics which is 
discussed in terms of the organic solid 6-ET2X. 

The spinless fermion Hamiltonian in the present study 
is given as, 



H 



-u 

(hj) 



h.c. + VijUifij^j. 



(1) 



Here, Cj and n j denote the annihilation and number op- 
erator of fermions. The triangular lattice has anisotropy 
in one of three directions as reflected in the transfer inte- 
grals, tij=t,t f , and the inter-site Coulomb interactions, 
Vij = V, V, between nearest neighbor sites whose corre- 
sponding bond directions are defined in Fig.l. The phe- 
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nomena realized in EHM in the temperature range of 
,2 

J < T<y are suitably described by the t-V model. 

We perform the exact diagonalization on N = 4 x 6 = 
24 cluster in eq.(l). Figure 1 shows the phase diagram 
on the plane of V/t and V /t at several fixed values of 
t'/t with periodic boundaries. The diagram is classified 
into three different regions, whose boundaries are esti- 
mated from the discontinuity in charge structural fac- 
tor, the gradient of energy by V and V ', and Drude-like 
component as we see shortly. The phases (I) and (III) 
are the two-fold periodic striped CO states whereas the 
phase (II) has a three fold periodicity and is denoted as 
"a pinball liquid" which is studied in Ref.[10]. The t' jt- 
dependence of the phase boundaries is very small. This 
phase diagram describes well the strong coupling physics 
at around V, V > 2t below which the phase transition 
begins to smear. 

We first present the structural factor, Ck = 
lEi m ((^-l/2)K-l/2))e^-^, in Fig. 2. Each 
phase is characterized by its own wave numbers, 
(I)(|-n, 7r) with integer n, (II) (7r, ±§tt), (0, =t~7r), and 
(III) (tt, 0) . The peak amplitudes vary discontinuously 
from one to another at the phase transitions, which 
confirm their first order character. We mention that 
the structural factors over the entire Brillouin zone be- 
side each characteristic ones are almost completely sup- 
pressed, which means that there is little possibility of 
the coexistence phase. The amplitudes in (I) and (II) 
are significantly smaller than the one in (III). As we see 
shortly, the disorder exists in the classical limit of phases 
(I) and (II), which still influences the states at t 7^ and 
suppresses the two-body correlation. 

Figure3(a) shows the ground state energy, Eq, which 
behaves almost linearly with V at the fixed value of 
V + V = 2Vb, while its gradient changes discontinu- 
ously at the phase transition. It is also seen in Fig. 3(b). 
Therefore the phase transitions are of first order. Note 
that the present phase diagram is almost free from finite 
size effect. Figure 3(a) reveals that Eq of the smaller 
cluster size, 4x4 and 4x3, reproduces within 0.3% the 
ones of 4 x 6 in (I), (III) and (II), (III), respectively, i.e. 
as long as the periodicity of each state is commensurate 
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Fig. 1. Ground state phase diagram of eq.(l) on the plane of 
V/t and V /t obtained numerically at N=2A. It is classified into 
three regions, (I) chain stripe, (II) pinball liquid and (III) vertical 
stripe. The phase boundary is given at t' /t=l, 0.5 and 0. Bro- 
ken lines with square symbols indicate the upper bound of the 
insulating state in (I) for t/t' = 0.5,1. Shaded region includes 
the possible location of 0-ET2AT. Representative configuration 
of the classical limit of the three phases are schematically shown 
together. The arrows with V, t and V ,t' at the top indicate the 
directions of interacting bonds in eq.(l). 



Fig. 3. Ground state energy of eq.(l), Eq : (a) at fixed V + V = 
2V = 18 with t = t' =1 under several cluster size. Broken lines 
give the classical binding energy per site, where E\,E2 are the 
energy gain per site due to t, t' . The 4x6 cluster shown together 
on the right, (b) ^-dependence of Eq at V=0. for several choices 
of t' = 0,0.5,1 (t = 1). Arrows indicate the phase transitions 
between (I) and (II). (c) Energy correction, (Eq —Ec)/N, due 
to t = t' and t = 2t' at V Q =2Q in three phases, V=8(I), 10(11), 
and 12(111), which are assigned E2, Ei, and E2, respectively. 
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Fig. 2. Amplitude of charge structural factors along V + V = 
2V =12 with characteristic wave numbers. Rectangle repre- 
sents the Brillouin zone and the circles are the location of wave 
numbers with peaks. 
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with the shape of the cluster. 

The difference between three phases can be clearly un- 
derstood by the strong coupling theory. 10 In the classi- 
cal limit t = t' = 0, the ratio of V and V determines 
the phase diagram. The system has a semimacroscopi- 
cally (~ 2^-fold) degenerate chain CO state and a non- 
degenerate vertical stripe CO state at V < V and V > 
V respectively, 10 and the phases (I) and (III) are their 
extensions. Their lowest order energy gain per site by t' 
(< V, V) is estimated as, E 2 = -t' 2 /V - t 2 /(2V' - V) 
(I), -2£ 2 /(3V-2V) (II), which is presented numerically 
in Fig. 3(a) as (Eq — Eq)/N. Since this effect is small 
at V, V ^> t, the stripes are quite robust even though 
the charge localization length is somewhat extended. 

When V = V the three sublattice state with a macro- 
scopic degeneracy joins the ground state besides these 



Fig. 4. (a) Three-body charge correlation function, P3C7), along 
the second chain(i = 7 — 12) at t = t' = 1 for several choices of V 
with V = 6, with i denotes the site index. Those in broken lines 
are the cases of striped states, (b) Fourth order process that lifts 
the degeneracy of the chain state. The cyclic permutation of three 
charges on the rectangle gives 64. The energy spectrum becomes 
equally spaced by e^Ny, with bottom and top as the horizontal 
and diagonal CO, respectively. The number kinks including A^- 
rectangles in bold lines determines the energy. 

two stripes in the classical limit. The quantum dynamics 
of t by first order lifts the degeneracy of part of the three 
sublattice states and the pinball liquid (II) appears. The 
estimated kinetic energy gain ((Eq — Ec) / 'N in Fig. 3(a)) 
is of first order (t-linear) which behaves as E\ ex cit+C2t 2 
(ci,C2 are constants) in contrast to E2 oc t 2 as shown 
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Fig. 5. Amplitude of energy, D a , under the twisted boundary 
condition in a = x,y, directions along (a)V=0 and (b)V=0. 
Solid lines show the exponential fit, D oi . Arrows 

indicate the phase transitions (those in (a) are the possible metal- 
insulator transition points). The corresponding fitting constants, 
C^ a , are given in (c) as functions of t' jt. 



in Fig. 3(c) regardless of the anisotropy in t 1 jt. Since 
E Q = V/2 + E 2 , (2V + V)/6 + E 1 , and V/2 + E 2 for 
(I)- (III), respectively, the phase (II) appears at around 
6(^1 - E 2 ) < V - V < -3(E 1 - E 2 ). 

To fully figure out the nature of the pinball liquid (II) 
out of disorder, we should consider the three-body cor- 
relation, Ps(j) = (ni(l — n 2 )rij) which is the popula- 
tion of the j-th site when the i-st and the 2-nd site is 
present and absent, respectively. Figure 4(a) shows that 
the ha + nB~2nc type of three- fold structure is present 
only in (II). This is in contrast to the tla ^ tlb ~ nc 
type ones based on mean field wave functions. 8 

The phase (I) is another ordered state out of disor- 
der. Its classical limit has a CO along every chain (y- 
direction). Regular stacking of the CO chains give the 
diagonal striped structure shown in Fig. 4(b). Starting 
from this we can introduce a kink one by one, and fi- 
nally we obtain the horizontal stripe with the full num- 
ber of kinks. All these structures have the same energy 
and are semi-macroscopically degenerate. We call this 
gapless and disordered state a chain CO. 10 The energy 
correction by t starts from the second order, but the 
degeneracy is not lifted until the fourth order; the de- 
gree of lifting per rectangular unit with three charges is 
e 4 = -2t A /(2V' - VfV as shown in Fig. 4(b). The en- 
ergy gain from the diagonal stripe is given as e^NyN^^ 
where A/kink denotes the number of kinks. Therefore the 
horizontal stripe becomes a unique ground state in the 
bulk limit, and the reduction of A/j^nk corresponds to 
the gapful domain wall excitation which requires energy 
larger than the single particle excitation (e 4 N y ^> V') ex- 
cept in the vicinity of the strong coupling limit. In fi- 
nite clusters, however, these kinked chain states easily 
mix with the horizontal stripe. The energy gain here is, 

Ny Ny 

e c h ~ N y t'^~ /V'^ - , which is far smaller than E\ and 
E 2 at V ^> t' so that its correction to the phase bound- 
ary in Fig. 1 is negligibly small. In this way, the degener- 
acy lifted ground state in the bulk limit should have only 
a single peak of Ck within the Brillouin zone at (0, tt). 
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Fig. 6. (a) Schematic ground state phase diagram on the planes 
of V/t and V/t at each t' It. The COI denote the insulating 
striped phases. Filled and open circles are the possible locations 
of the Rb and Cs-salts of 0-ET2X. (b) Schematic phase diagram 
at finite temperature. The colored regions denote the quantum 
phases, (I)- (III), which are replaced by the classical disordered 
phase at higher temperatures. The crossover is possibly realized 
in the Cs-salt at T cr as indicated in arrows. 



Next, we twist the boundary condition in x and y di- 
rections by the phase factor of e 1 ^ (a = x,y) : where 
^=0 and 7r correspond to the periodic and anti-periodic 
boundary conditions, respectively. Then, D a = EQ(ip a = 
7r) — EQ(ip a =0), is approximately proportional to the 
Drude weight of the finite system size. The discontinu- 
ity of D a detects the first order phase transitions among 
(I) (II) and (III). The pinball liquid is possibly a metal 
in the bulk limit with D a ^ as previously reported in 
Ref.[10]. 10 In contrast, the coherence of the gapful local- 
ized commensurate stripes is expected to decay exponen- 
tially with distance I as, ex e _z /^, where £ denotes the co- 
herence length. Then we except, £ a = C^ a /A c with C^ a 
as V, V'-independent but ^-dependent values and A c 
the value of charge gap in the classical limit, i.e. V in (II) 
and 2V - V in (III). By anticipating D a ex e ~ NaC ^ /Ac 
(N a is the system length), the fitting is well performed 
as the solid lines in Figs. 5(a) and 5(b) whose Cgjs are 
plotted against t' jt in Fig. 5(c). The fitting breaks down 
simultaneously at the transition from (III) to (II), which 
indicates that phase (III) is a localized (£ a < N a ) insula- 
tor. In this phase both D a and C^ a are t'/t-independent 
which is also the case with Eq in Fig. 3(b). In contrast, 
Eq of (I) in Fig. 3(d) together with D a actually has 
large ^/t-dependence. There is a possible phase bound- 
ary within phase (I) estimated by the discontinuity in 
dD a /dV which is given in Fig. 1 with broken lines. The 
strong coupling large V'-region is an insulator in its bulk 
limit. However, in the rather weak coupling region, £ a de- 
velops and exceeds which corresponds to the above 
phase boundary. Notice that it depends much out' jt. At 
present, we cannot discriminate whether this new inter- 
mediate region in (I) is gapful or not in the bulk case, nor 
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whether there is a finite size effect in this boundary. Nev- 
ertheless, the estimated boundary gives an upper bound 
of the metal insulator transition line. 

The previous theoretical attempts to clarify the phase 
diagram suffers from problems. Those based on mean- 
field wave function 8 failed to describe fully the extremely 
correlated chain(horizontal) and pinball states. Their 
phase diagram has ambiguous "competition" region be- 
tween states of different symmetry breaking due to large 
size-dependence, which is an artifact of adopting a plane 
wave based wave function on such correlated and highly 
localized system. Exact diagonalization on EHM in 4x4 
finite clusters overestimates their "frustration induced 
metal" as well as underestimates the pinball liquid, 9 as 
clarified explicitly in Figs. 1 and 3(a). Our phase diagram 
presents a clear view to these problems and describes the 
intrinsic nature of U = oo EHM at quarter-filling since the 
basic character of three phases are free from statistics. 10 

The ground state phase diagram is summarized in 
Fig. 6(a). The anisotropy of t' /t only affects the coher- 
ence of phase (I) in a relatively weakly coupled region. 
We show schematically in Fig. 6(b) the finite tempera- 
ture effect, fe^T^O, where the classical disordered state 
becomes prosperous as a result of the macroscopic en- 
tropy. The free energy of the classical three-sublattice 
state is, F 3sub ~ 2y + y ' N-k B T ln(0.5 x 2 N ), while those 
of the nondegenerate quantum states are F = Eq ~ 
E c + SEqN with SE Q ~E 2 + e 4 , E u and E 2 for (I), (II) 
and (HI), respectively. The crossover lines from quan- 
tum to classical states with increasing ksT are esti- 
mated by the crossing points of these free energies as, 
T CT ~ -E 2 ,\E 1 \ and - E 2 , respectively. The 

chain type of disorder does not appear in the bulk case 
since the entropy gain remains of order ~ 0(\ /f N In y/N). 

Now we examine the relevance of the present results 
to the CO materials, 0-ET 2 X. The nature of this fam- 
ily at low temperatures is classified on the experimen- 
tal phase diagram along the dihedral angle, when 
(j) > 110°, the system becomes a CO insulator accom- 
panied by the structural transition at Tqo ~ 200K. At 
(j) < 110° it is a bad metal down to low temperatures 
with the resistivity minimum at T p ~ 20-70K. Recent 
X-ray analysis reports the existence of a short range 
diffuse spot of two-fold and three-fold type in the Cs- 
salt (X=CsM(SCN) 4 , M=Co,Zn, < 110°) at - T p 13 
and also in Rb-salts(0 ~ 110°) above Tco, 14 which have 
(\t'\/t,V'/V) ~ (0,0.9) and (0.5,0.85), respectively. 15 

It is obvious that there is a characteristic crossover 
temperature of order T p at <p < 110° which increases 
with (j). Actually, at around V > V, we can estimate the 
crossover temperature as T cr ~ v ~ v ~ 0.1V. Since we 
consider V ~ Tqo, the experimental crossover at ~ T p 
is reasonably understood. An observed bad metallic na- 
ture with almost T-independent resistivity at T > T p 12 
supports the existence of disorder. Also the up going of 
resistivity at T < T p is consistent with the insulating 
character of the possible horizontal stripe. 12 The coex- 
istence of states with different periodicity in the ground 
state is ruled out in the present study. However, at finite 
temperature, the disordered state which includes both 
the three-sublattice and the striped local patterns might 



appear around the crossover temperature, which were 
excluded in Fig. 6(b) for simplicity. The study on this 
complicated matter is still underway. 

Somehow, the ground state of this family is placed 
amidst the pinball liquid state as shown in Fig. 1, which is 
inconsistent with the above experimental findings. In the 
present stage we did not consider the effect of phonons. 
The three-fold state is found to be quite easily replaced 
by the horizontal stripe long range order(LRO) under 
that effect, 17 which might be valid in our case as well. 

It is frequently argued that (j) can be regarded as V/t. 
Contrarily, one of the authors successfully reclassified the 
experimental phase diagram by \t f \/t, on the basis of 
EHM 7 which is recently pursued experimentally 16 and 
also in the weak coupling theory. 8 In our phase diagram, 
however, \t'\/t only rules the coherence length of the sys- 
tem at the relatively weak coupling region. To figure this 
out, we should examine the extra effect such as electron- 
phonon coupling in the disordered region. 

In conclusion, the strong coupling phase diagram of 
the charges on the anisotropic triangular lattice is pre- 
sented. A metallic pinball liquid phase appears in be- 
tween an insulating vertical and the horizontal stripe 
CO's. The first order transitions take place in the phase 
boundaries which are almost solely determined by the 
geometry of the inter-site Coulomb interactions. The 
present phase diagram provides two different types of 
order-from-disorder due to quantum dynamics; the ex- 
tremely correlated pinball liquid out of disordered three 
sublattice classical state, and the horizontal state which 
has a domain- wall like excitation to the disordered chain 
striped classical states. Both of these states provide good 
test cases for the disorder-embedded phenomena. The 
disorder once gave path to quantum states recovers at fi- 
nite temperature with the aid of its macroscopic entropy. 

We thank R. Kondo, M. Watanabe, I. Terasaki, and 
R. Chiba for many useful experimental informations. 
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